clc;%清理命令行窗口
data=readmatrix('RECposition.txt');%读取excel数据，同目录下
ECEF_X=-2612403.257446;
ECEF_Y=4740776.195834;
ECEF_Z=3362563.549283;
figure(1);% Longitude/latitude/height trajectory lines
plot3(data(:,5),data(:,6),data(:,7),'k:.');
title('BLH trajectory lines');xlabel("longtitude/°");ylabel("latitude/°");zlabel("high/m");
hold on;
figure(2);% Longitude/latitude trajectory lines
plot(data(:,5),data(:,6),'k:.');
title('Latitude/longitude line');xlabel("longtitude/°");ylabel("latitude/°");
hold on;
figure(3);% Longitude/latitude height Curve
subplot(1,3,1);
plot(data(:,5),'k-');
xlabel("历元数/×0.1s");ylabel("longtitude/°");
subplot(1,3,2);
plot(data(:,6),'k-');
xlabel("历元数/×0.1s");ylabel("latitude/°");
subplot(1,3,3);
plot(data(:,7),'k-');
xlabel("历元数/×0.1s");ylabel("high/m");
hold on;
figure(4);% ECEF Coordinates X Y Z Curve
subplot(1,3,1);
plot(data(:,2),'k-');
xlabel("历元数");ylabel("ECEF-X/m");
subplot(1,3,2);
plot(data(:,3),'k-');
xlabel("历元数");ylabel("ECEF-Y/m");
subplot(1,3,3);
plot(data(:,4),'k-');
xlabel("历元数");ylabel("ECEF-Z/m");
hold on;
figure(5);% clock bias
subplot(1,3,1);
plot(data(:,8),'k');
xlabel("历元数");ylabel("reciver CLOCK BIAS/m");
subplot(1,3,2);
plot(data(:,9),'k');
xlabel("历元数");ylabel("Galileo-GPS/m");
subplot(1,3,3);
plot(data(:,10),'k');
xlabel("历元数");ylabel("BDS-GPS/m");
hold on;
figure(6);% Dilution of Precision(DOP)
subplot(1,4,1);
plot(data(:,11),'k-');
xlim([0 540081]);
ylim([0 3]);
xlabel("epoch");ylabel("GDOP(Geometric Dilution of Precision)");
subplot(1,4,2);
plot(data(:,12),'k-');
xlabel("epoch");ylabel("PDOP(Position Dilution of Precision)");
subplot(1,4,3);
plot(data(:,13),'k-');
xlabel("epoch");ylabel("HDOP(horizontal dilution of precision)");
subplot(1,4,4);
plot(data(:,14),'k-');
xlabel("epoch");ylabel("VDOP(vertical dilution of precision)");
figure(7);% ECEF Coordinates X Y Z Error
subplot(1,3,1);
plot((data(:,2)-ECEF_X)/10,'k-');
xlabel("time/h");ylabel("ECEF-X error/m");
xlim([0 540081]);
subplot(1,3,2);
plot((data(:,3)-ECEF_Y)/10,'k-');
xlabel("time/h");ylabel("ECEF-Y error/m");
xlim([0 540081]);
subplot(1,3,3);
plot((data(:,4)-ECEF_Z)/10,'k-');
xlabel("time/h");ylabel("ECEF-Z error/m");
xlim([0 540081]);
hold on;